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Abstract 

Planar Cell Polarity (PCP) is an evolutionarily conserved characteristic of animal tissues marked by coordinated polarization 
of cells or structures in the plane of a tissue. In insect wing epithelium, for instance, PCP is characterized by en masse 
orientation of hairs orthogonal to its apical-basal axis and pointing along the proximal-distal axis of the organ. Directional 
cue for PCP has been proposed to be generated by complex sets of interactions amongst three proteins - Fat (Ft), Dachsous 
(Ds) and Four-jointed (Fj). Ft and Ds are two atypical cadherins, which are phosphorylated by Fj, a Golgi kinase. Ft and Ds 
from adjacent cells bind heterophilically via their tandem cadherin repeats, and their binding affinities are regulated by Fj. 
Further, in the wing epithelium, sub-cellular levels of Ft-Ds heterodimers are seen to be elevated at the distal edges of 
individual cells, prefiguring their PCP. Mechanisms generating this sub-cellular asymmetry of Ft-Ds heterodimer in proximal 
and distal edges of cells, however, have not been resolved yet. Using a mathematical modeling approach, here we provide a 
framework for generation of this sub-cellular asymmetry of Ft-Ds heterodimer. First, we explain how the known interactions 
within Ft-Ds-Fj system translate into sub-cellular asymmetry of Ft-Ds heterodimer. Second, we show that this asymmetric 
localization of Ft-Ds heterodimer is lost when tissue-level gradient of Fj is flattened, or when phosphorylation of Ft by Fj is 
abolished, but not when tissue-level gradient of Ds is flattened or when phosphorylation of Ds is abrogated. Finally, we 
show that distal enrichment of Ds also amplifies Ft-Ds asymmetry. These observations reveal that gradient of Fj expression, 
phosphorylation of Ft by Fj and sub-cellular distal accumulation of Ds are three critical elements required for generating 
sub-cellular asymmetry of Ft-Ds heterodimer. Our model integrates the known experimental data and presents testable 
predictions for future studies. 
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Introduction 

Planar Cell Polarity (PCP) represents coordinated orientation of 
cells or structures in an axis within the tissue plane. Most animal 
organs display PCP, which when perturbed lead to developmental 
anomalies such as open neural tubes, polycystic kidneys, and even 
can cause deafness due to loss of planar polarization of the cilia in 
the inner ear [1]. PCP has been extensively studied in Drosophila 
wing, eye, and in the adult abdomen [2]. Drosophila wing, in 
particular, is an ideal model system to study PCP, which is 
mirrored by projection of actin-rich trichomes (hairs) from the 
distal edges of its individual cells (Figure ID) [3,4]. 

Three evolutionarily conserved proteins, namely, Fat (Ft), 
Dachsous (Ds) and Four-jointed (Fj), the so called Ft-Ds-Fj system 
or global regulators of PCP, generate a tissue-level directional cue 
thereby linking the axis of polarization of individual cells with that 
of the tissue/organ: for instance, projection of the wing hairs along 
the proximal-distal axis in Drosophila wing [5]. A second class of 
proteins consisting of Frizzled (Fz), Flamingo (Fmi), Van Gogh 



(Vang), Prickled (Pk) and Dishevelled (Dsh), referred to as core 
PCP proteins, propagate PCP signals from cell-to-cell [2,6]. 
Finally, several down-stream effectors, including cytoskeletal 
proteins, respond to these signals from global regulators and core 
proteins to establish PCP [7]. Upon loss of global signal, as in ft, ds, 
Jj mutants, PCP could be seen to be coordinated amongst 
neighboring cells; however, the direction of their polarity is often 
not aligned with the tissue/organ axis thereby giving rise to 
swirling patterns (Figure ID) [8- 10]. Loss of activities of the core 
proteins, in contrast, abolish PCP altogether in individual cells; for 
instance, wing hairs in fz mutant may grow out of the centers of 
the cells, instead of their distal edges [6] . 

Ft, Ds and Fj display intricate interactions with each other. Ft 
and Ds, for instance, are large atypical cadherins that bind 
heterophilically between adjacent cells via their tandem cadherin 
repeats [5]. Fj, a golgi kinase, phosphorylates both Ft and Ds 
[11,12], which in turn increases affinity of Ft for Ds and decreases 
that of Ds for Ft [13,14]. Also, expression levels of Ds and Fj vary 
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across the tissue; in the wing primordium, for instance, these form 
opposing gradients where Fj expression is higher in its distal than 
that of its proximal domain, and vice-versa for Ds (Figure 1A) 
[10,15,16]. These opposing gradients of Fj and Ds expressions 
have been seen to be critical for establishing PCP [9,17]. Further, 
Ft-Ds heterodimer also displays a sub-cellular distal enrichment, 
thereby, resulting in its asymmetric localization (polarization) in 
the proximal and distal membrane of the cells preceding their 
overt display of PCP through asymmetric localization of the core 
proteins such as Fz [18,19]. Ft-Ds heterodimer asymmetry thus 
could be the earliest and coarse signals for cell polarization that are 
eventually strengthened by sub-cellular asymmetry of the core 
PCP proteins and their signaling. A critical question that has not 
been answered yet is how the known interactions within these 
three elements of the Ft-Ds-Fj system translate into a sub-cellular 
asymmetry of the Ft-Ds heterodimer. 

A mathematical modeling approach can offer insights into how 
a set of biochemical interactions at a local (sub-cellular) scale could 
impact a global pattern such as the coordinated PCP in an entire 
organ like the insect wing. Several mathematical models have been 
developed to understand various aspects of PCP [20] . The first 
integrated modeling and experimental study in PCP [21] offered a 
model for intercellular feedback regulations amongst core proteins 
to explain 'domineering effect' [22] - a phenomenon wherein 
perturbation in PCP in a group of cells lacking (mutant), for 
instance, Fz protein activity was relayed to wild type neighboring 
cells. This study proposed that domineering effect can simply 
emerge from the dynamics of the known intercellular and 
intracellular interactions among the core proteins (Fz, Vang, Pk 
and Dsh). This model further suggested that a directional cue 
provided by Ft-Ds-Fj system is crucial for explaining the observed 
PCP patterns [21]; however, it did not elucidate how the Ft-Ds-Fj 
system as such generates this cue. Further, attempts at mathemat- 
ical modeling of PCP captured the mechanism of action of the 
core proteins, and strengthened the idea that a weak tissue-level 
directional cue is important for obtaining the experimentally 
observed PCP patterns [23— 30]. A recent phenomenological model 
proposes how Ft-Ds heterodimer formation in a group of cells in 
the tissue can induce collective polarization over the entire tissue; 
however, this model did not factor the crucial question of the role 
played by tissue-level expression gradients of Ds and Fj for the 
regulation of PCP [31]. Since none of these models factored all 
activities of the members of Ft-Ds-Fj system, our current 
understanding of their contribution in PCP regulation remains 
far from complete. 

Here, using a mathematical modeling approach, we have 
simulated the experimentally identified interactions amongst the 
members of Ft-Ds-Fj system to reveal their impact on sub-cellular 
asymmetry of Ft-Ds heterodimer. We have developed an Ordinary 
Differential Equation (ODE)-based model by taking into consid- 
eration the known interactions in Ft-Ds-Fj system. These elements 
are: (a) Ft-Ds heterophilic binding, (b) opposing tissue-level 
gradients of Ds and Fj and, finally, (c) phosphorylation of Ft and 
Ds by Fj kinase. We show that Fj gradient, phosphorylation of Ft 
by Fj, and, sub-cellular distal accumulations of Ds are three 
essentials for establishing sub-cellular asymmetry of Ft-Ds hetero- 
dimer in individual cells. 

Results and Discussion 

Model Formulation 

We assumed that cells in Drosophila wing are arranged along 
rows that are oriented in a single proximal-to-distal axis. We 
modeled one such row of cells for deciphering the sub-cellular 
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Figure 1. Fat (Ft)-Dachsous (Ds)-Four-jointed (Fj) system 
regulates PCP in Drosophila wing. (A) Schematic of Drosophila wing 
showing opposing expression gradients of Ds (green) and Fj (grey) 
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while expression of Ft (not shown) remains uniform. Arrow marks the 
proximal-distal (P/D) axis of the wing. (B) Arrangement of cells (boxes) 
seen in the contexts of proximo distally increasing gradient of Fj 
expression (grey) as well its activity; the latter represented by arrows 
and hammers symbolizing its activating and repressive roles, respec- 
tively, while thicknesses of these lines indicate relative strengths of 
these activities. (C) Different concentrations of the Ft-Ds heterodimer 
(represented by thickness of red line) at the proximal and distal edges 
of the cells provide distinct intensity plot. (C) A line joining the peaks of 
these two levels of Ft-Ds heterodimer at the proximal and distal edges 
would be proximally or distally inclined depending on the concentra- 
tions of the Ft-Ds heterodimer at these edges. This line would, however, 
be horizontally placed when concentrations of the Ft-Ds heterodimer 
are equal at the proximal and distal edges. (C") Simplified visual 
representations of the concentrations of Ft-Ds heterodimer at the 
proximal and distal edges of the epithelial cells seen in (C). (D) 
Orientation of wing hairs in wild type and in dachsous (ds) mutant adult 
wings. 

doi:1 0.1 371 /journal.pone.0097641 .g001 

asymmetry of Ft-Ds heterodimer (Figure IB). For the ease of visual 
representation, we have denoted asymmetry in sub-cellular 
enrichment of Ft-Ds heterodimer in proximal and distal edges of 
individual cells by trapezoid shapes where the slope, either distally 
or proximally inclined, mirrors its enrichment at the distal or 
proximal cell membrane, respectively. However, when Ft-Ds 
heterodimers are equally enriched in proximal and distal edges, 
their symmetric distribution is represented in a rectangular shape 
(Figure 1C, C and C"). We have assumed the physical gradient of 
Fj to mean a gradient for its kinase activity; thus, in our 
simulations, the Fj kinase activity gradient is synonymous with its 
actual spatial gradient of expression (Figure 1A). Further, due to a 
distally enriched graded expression of Fj in Drosophila wing, its 
kinase activity too is shown to be graded along the proximal-distal 
axis of the wing (Figure IB). 

Since Ds and Fj are present in opposing gradients in the wing, 
total concentration of Ds (represented by Ps° for /* cell) decreases 
from proximal to distal end of the wing. On the other hand, total 
concentration of Ft (represented by Ff) is taken to be same in all 
cells. We assume that during establishment of PCP, total 
concentrations of Ft and Ds (sum of both phosphorylated forms 
of Ft and Ds) remain uniform in each cell. However, due to the 
gradient of Fj kinase activity, levels of phosphorylated Ft and Ds 
turn out to be different for each cell. Thus, the only variables 
considered during establishment of PCP are phosphorylated and 
unphosphorylated states of Ft and Ds. Fj phosphorylates both Ds 
and Ft and their rates of phosphorylation are designated ay and 
a.j, respectively, while their dephosphorylation rates are designat- 
ed as jSy and P d , respectively. The effect of Fj on phosphorylation 
of Ft and Ds in i h cell (Figure IB) is considered to follow a 
monotonically increasing Hill's functionG(iy'i)) which is given by 



G(Fj,) = 



Fj n + K" 



(1) 



dFf t 
dt 



-OfFifGiFfii-pfFf, 



(2) 



dPs p 
dt 



= a d Ps u i G{Fj i )-p d Ds p i 



(3) 



where Fji,Ff f ,Ds^ , Ft" and Ps" represent the concentrations of Fj, 
phosphorylated Ft, phosphorylated Ds, unphosphorylated Ft, and 
unphosphorylated Ds, respectively in the 2 th cell. 

We performed the analysis with randomly initialized values of 
concentrations of phosphorylated Ft and Ds (Figure 2A-H). In the 
SI (File SI), we show that the levels of phosphorylated Ft (Ff) and 
Ds (Ps p ) attain stable steady state values over time irrespective of 
the initial values of Ft p and Ps p . Thus, our results hold good for 
any set of initial values of Ff and Ps p in the row of cells 
considered. 

In our analysis, we normalized the concentrations of phosphor- 
ylated and unphosphorylated Ft by/ 7 ? 0 (the total concentration of 
Ft in each cell). Also, we normalized the levels of Ds and Fj with 
their respective total concentrations in the most distal cell of the 
row. We can calculate the concentration of unphosphorylated Ds 
in the i h cell as 



Ds" = Ps°-P^. 



(4) 



Since phosphorylated Ds (Ps p ) has decreased affinity to bind to 
Ft, while phosphorylated Ft (Ft p ) has increased affinity to bind to 
Ds [13,14], we assumed that the interaction between phosphor- 
ylated Ft and unphosphorylated Ds (Ft p — Ps") is much greater 
than those between other forms of Ds and Ft (i.e. Ft p —Ps p , 
Ff — Ds" and Ft"-Ps p ). Thus, Ft p -Ps" is the most abundant 
form of Ft-Ds heterodimer. Therefore, hereafter, by Ft-Ds 
heterodimer, we refer to the binding of phosphorylated Ft,Ft p , 
of one cell with unphosphorylated Ds, Ps", of its adjacent cell. 
Thus, two different binding events contribute to the formation of 
Ft-Ds heterodimer at the edge between two given cells A and B - 
(a) Ff of A binds with Ps" of B, and (b) vice-versa i.e., Ff of B 
binds with Ps" of A. Concentration of this Ft-Ds heterodimer is 
referred to as C. If the rates of formation and dissociation of the 
heterodimers are a c and /? Ci respectively, then the rates of 
heterodimer formation at the two edges of the i" 1 cell are given 
by equations (5) and (6) 



dC 
dt 



('-1.0 



= a c {Ff i Ps" l _ l +Ps" i Ff i _ x )-P c C\, l 



(i-l,/)> 



(5) 



where n and K are two model parameters. 

Since Fj is expressed in a tissue-level gradient, its kinase activity- 
is considered to increase linearly as we move along the proximal- 
distal axis of the wing. Hill's function thus ensures that effect of Fj 
on phosphorylation of Ft and Ds is higher in the distal end as 
compared to that in the proximal. This difference in the strength 
of regulation is shown by the increasing thickness of arrows and 
hammers in Figure IB as we move distally along the row of cells. 
Therefore, we can write 



dC 
dt 



(M+l) 



= * c {Ff l Ps1 +l+ Pf i Ff i+x )-P c C\ u+l) . 



(6) 



Equation (5) represents the heterodimer formation at the 
proximal edge of i cell, or the edge shared between i h cell and 
(i—lf l cell, whereas Equation (6) represents the heterodimer 
formation at the distal edge of i h cell, or the edge shared between 
i h cell and (i+l)' k cell. Therefore, the rate of evolution of sub- 
cellular asymmetry in the heterodimer formation between the two 
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Cell number 

Figure 2. Evolution of phosphorylated Ft [Ff) and phosphorylated Ds (Df). (A-D) Evolution of sub-cellular concentrations of 
phosphorylated Ft (Ff) and that of (E-H) phosphorylated Ds (Ds p ) in cells over time. (A) and (E) show the randomized initial concentrations of Ft 1 ' 
and Ds p , respectively. 
doi:1 0.1 371 /journal.pone.0097641 .g002 



edges of the i' 1 cell is given by equation (7). 



dC 



dC 

Of, 



(i 1,0 



Ft, (D$ + , - Ds»_ , ) + Ds" {Ff i+ 
~ Pc\C\(i,i+\) ~ Q(i-l,i)J 



(7) 



Furthermore, we assumed that the heterodimer formation does 
not influence the levels of available Ft and Ds monomers 
(Ft" ,Ds" ,Ft p ,Ds p ) in the cell. Also, we did not consider any 
feedback on monomer concentration due to heterodimer forma- 
tion, since heterodimer formation and phosphorylation of unpho- 
sphorylated monomers (Ft" ,Ds") happens in different cellular 
compartments - the former at cell edges, and the latter at the 
Golgi by the Fj kinase. 

Parameter Values and Initial Conditions 

Representative model parameter values for the simulations are 
listed in Table 1. The sensitivity analysis (see Figures SI, S2 in File 
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SI) showed that the conclusions drawn from the model do not 
change for a large range of parameter values. 

Ft-Ds-Fj Interactions Culminate in Sub-cellular 
Asymmetric Enrichment of Ft-Ds Heterodimer at the 
Distal Edges 

The model developed was numerically simulated for 20 cells. 
Simulations for phosphorylated Ft,Ft p (Figure 2 A-D), and 
phosphorylated Ds,Z) 1 ? / '(Figure 2 E-H), show that their levels in 
a cell achieve a steady state over time. The expression for time 
evolution of Ftf and Ds p is given in SI (File SI). We also simulated 
the model for many different sets of parameters, and found quite 
similar results for the set of parameters as present in Table 1. 
Variations in the parameters, namely, n, K, ay^/^OQ/^displayed 
only a modest effect on our results for a very large range (Figures 
SI, S2 in File SI). We also varied the relative strength of the kinase 
activity of Fj for its substrates, Ft and Ds, by changing the 
parameters of the two Hill's functions in equations (2) and (3), and 
found our results not to be overtly sensitive to these changes too 
(Figure S2). 

Concentration of Ft-Ds heterodimer (or the Ft p — Ds" hetero- 
dimer) across the tissue evolved with time, starting with random 
uncoordinated values to a coordinated steady state at the cell edges 
(Figure 3 A-H). In all panels of Figure 3, each trapezoidal column 
represents levels of Ft-Ds heterodimer on the two edges of a single 
cell, where the apparent slope of the trapezoid, as explained in 
Figure 1 C, C and C", mirrors disparity in the concentrations of 
Ft-Ds heterodimer at proximal and distal cell edges, revealing their 
asymmetric sub-cellular enrichments. Thus, in a steady state, 
concentration of Ft-Ds heterodimer for each cell in wild type 
scenario at its distal edge is more than that on its proximal edge 
(Figure 3H). Also, since distal edge of a cell abuts the proximal 
edge of its distal neighbor, value of Ft-Ds heterodimer concentra- 
tion at the distal edge for one cell equals that of the proximal edge 
of its neighbor (Figure 3H). 

These simulations thus show that interactions within the Ft-Ds- 
Fj system culminate in a sub-cellular distal enrichment of Ft-Ds 
heterodimer (Figure 3H). This is corroborated by recent exper- 



PLOS ONE | www.plosone.org 



4 



May 2014 | Volume 9 | Issue 5 | e97641 



Sub-Cellular Asymmetry of Ft-Ds Heterodimer for Planar Cell Polarity 



imental observations in the wing that display sub-cellular 
asymmetric localizations of the Ft-Ds heterodimer [18,19]. In 
these studies, Ft and Ds were seen to localize asymmetrically on 
the two edges of the cell, with low Ft and high Ds on the distal 
edge, and high Ft and low Ds on the proximal [18,19]. 

The model presented here appears to be robust in generating 
consistent polarity patterns over a large range of parameter values, 
and variations in initial conditions of randomized Ft and Ds 
activity (see SI File SI). Besides a sub-cellular graded pattern of Ft- 
Ds heterodimer, the model also predicts a distally elevated 
endogenous tissue-level graded pattern of Ft-Ds heterodimer 
formation (Figure 3H), which has not yet been observed 
experimentally. This spatial gradient in Ft-Ds heterodimer 
concentration is likely to have domain-specific impact on PCP of 
the organ. Uniformly high levels of Ft or Ds, on the other hand, 
would alter this graded pattern of Ft-Ds heterodimer, which in 
turn could preferentially affect PCP in the proximal, but not the 
distal wing - endogenous concentration Ft-Ds heterodimer being 
low in the former. Previous experimental results appear to 
corroborate this fallout of our simulation. Thus, when Ft or Ds 
was overexpressed in their respective naive contexts: that is, in ft or 
ds mutant background, PCP in the distal but not in proximal 
domain was rescued [16,32]. 

Also, since the Ft-Ds heterodimer represents phosphorylated Ft 
(Ft p ) bound with unphosphorylated Ds (Ds") (see Model formu- 
lation), a sub-cellular asymmetry (distal enrichment) of the 
heterodimer would also mean distally enhanced phosphorylated 
Ft (Ft p ). It may be noted that previous speculation [9,16] about 
asymmetric Ft activity (which is synonymous with its phosphor- 
ylated state) has not been validated yet. Immunochemical 
detection using antibody that discriminates between phosphory- 
lated and unphosphorylated Ft could provide an assay for distally 
enhanced (asymmetric) Ft activity. 



Fj Gradient, but not Ds Gradient, Plays an Essential Role 
in Attaining Sub-cellular Asymmetry of Ft-Ds 
Heterodimer 

To test if the distally elevated gradient of Fj at the tissue level in 
the developing wing [11,16] is critical for attaining sub-cellular 
asymmetric localization of the Ft-Ds heterodimer, we examined 
the consequences of flattening Fj gradient in silico. Indeed, when 
tissue-level expression gradient (or activity gradient; since we 
assume the expression and kinase gradients of Fj to be 
synonymous) of Fj was abolished by its uniform overexpression, 
it led to a loss of sub-cellular asymmetry in Ft-Ds heterodimer 
localization (Figure 4B). We expect this loss of sub-cellular 
asymmetry of Ft-Ds heterodimer to cause the loss of sub-cellular 
asymmetry of developing wing cells. Our findings corroborated 
well with the experimental observations wherein overexpression of 
Fj (leveled Fj gradient) in the wing leads to the loss of membrane 
localization of the an unconventional myosin Dachs (D) [33]. It 
may be noted that D is a downstream component of the Ft-Ds-Fj 
system, which co-localizes with Ds on the distal edges of the cells 
[33], and thus presents a readout of sub-cellular asymmetry and 
cell polarization [18,19,34]. 

Like Fj, Ds too is expressed in a tissue-level gradient, but with a 
proximal elevation [10,16]. We examined the consequences of 
flattening this Ds gradient by its overexpression, and found that a 
flattened Ds gradient does not disturb sub-cellular asymmetry of 
Ft-Ds heterodimer, as long as the Fj gradient is intact (Figure 4C). 
This would imply that when Fj displays tissue-level gradient, 
expression of Ds alone - irrespective of its expression being 
uniform or in gradient at a tissue level - suffice for PCP. We then 
predicted that uniform overexpression of Ds would rescue polarity 
defects, for instance, in ds mutant wings since their endogenous Fj 
gradient remains unaffected; this was indeed observed experimen- 
tally [16]. On the other hand, a flattened Fj gradient disturbs sub- 
cellular asymmetry of the heterodimer, even when Ds gradient is 
intact (Figure 4B). Further, when both Ds and Fj gradients are 
flattened, the pattern observed (Figure 4D) appears akin to that 
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Figure 3. Evolution of asymmetric localization of the Ft-Ds heterodimer in cells over time. (A-H) Evolution of the asymmetric levels of Ft- 
Ds heterodimer in cells. (A) shows the initial uncoordinated levels of Ft-Ds heterodimer while that shown in (H) represent the final steady state. 
doi:1 0.1 371 /journal.pone.0097641 .g003 
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Figure 4. Fj gradient, but not that of Ds, is essential for generating the asymmetric localization of the Ft-Ds heterodimer. (A) 

Asymmetric localization of Ft-Ds heterodimer (distally heightened top of trapezoid) is seen when both- Fj and Ds are expressed in gradient. Data 
shown here is the same as that shown in Figure 3H. (B) Uniform over-expression of Fj results in a loss of asymmetry of Ft-Ds heterodimer (flattened 
top) in all cells, although Ds expression is maintained in a gradient. (C) Conversely, uniform overexpression of Ds does not affect Ft-Ds heterodimer 
asymmetry when Fj expression remains in a gradient. (D) When both Fj and Ds are uniformly overexpressed, asymmetry of Ft-Ds heterodimer is lost. 
doi:1 0.1 371 /journal.pone.0097641 .g004 



caused by flattening of Fj (Figure 4 B) but not the one presented by 
flattened Ds gradient (Figure 4C). Therefore, Fj gradient, but not 
that of Ds, is essential for attaining sub-cellular asymmetry of Ft- 
Ds heterodimer. In other words, Ds gradient by itself (without the 
Fj gradient) is not enough to give correct patterning. Together, 
these findings suggest that while Ds is required for generation of 
Ft-Ds asymmetry and cell polarization, its tissue-level gradient per 
se may not be critical for PCP. 

Phosphorylation of Ft by Fj is Essential in Achieving Sub- 
cellular Asymmetry of Ft-Ds 

Fj kinase phosphorylates both Ft and Ds and regulates binding 
of Ft to Ds in two ways - by inhibiting binding of Ds with Ft on 
one hand and by promoting binding of Ft with Ds on the other 
(Figure 5A) [13,14]. This opposing action of Fj on Ds and Ft has 
been proposed to establish sub-cellular polarity with high fidelity 
[6,14]. We asked whether the phosphorylation of both Ft and Ds 
by Fj kinase is necessary for producing sub-cellular Ft-Ds 
asymmetry. In silico mutation, which retained Fj kinase activity 
(phosphorylation) against Ds but lacked that against Ft (Fj' Ft ' 
mutant) disrupted PCP across the tissue (Figure 5B). The converse, 
meaning intact Fj kinase activity against Ft but loss of that against 
Ds (Fj (Ds) mutant), however, did not disrupt PCP (Figure 5C). 
Thus, while phosphorylation of Ft by Fj kinase appears to be an 
essential for achieving Ft-Ds sub-cellular asymmetry, redundancy 



of phosphorylation of Ds by Fj as seen in our simulation may 
confer robustness to polarization process across the tissue. 

Sub-cellular Polarity of Unphosphorylated Ds (Ds u ) 
Amplifies Sub-cellular Asymmetry of Ft-Ds Heterodimer 

Distal edge of a cell abuts proximal edge of its distal neighbor 
or, in other words, the edge shared between i"' cell and (i+lf cell 
is distal for i' k cell but proximal for (i+lf cell. At this edge, Ft-Ds 
heterodimer can be formed in two different orientations-Ds" 
(unphosphorylated Ds) of the i h cell binds with Ftf (phosphory- 
lated Ft) of the (i+lf cell (denoted by Ff t Ds" [+X ), ori^of the i' h 
cell binds to Ds"o{ (i+lf cell (denoted by D^Ff j+x ), as shown in 
Equation (6). Similarly, heterodimers formed at proximal edge of 
the i h cell, or the edge shared between t h and (i— if cell, can be 
formed in two possible orientations - Ds"oi the i ih cell binds with 
FtP of the (i-lf ceU (denoted by DsJF^^) or Ff of the t* ceU 
binds to Ds" of the (i-lf cell (denoted by F^D^_{), as shown in 
Equation (5). 

We calculate the difference in heterodimer concentration 
between distal and proximal edges of the i ih cell at steady state, 
denoted by AC* 9 and given by the equation given below (For 
details, see equation (S4) in File Sl)- 
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10 
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10 
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20 



Figure 5. Phosphorylation of Ft by Fj kinase is necessary for Ft-Ds asymmetry. (A) Wild type Fj activity displays characteristic distal 
enrichment of Ft-Ds heterodimers. Data shown here is same as that shown in Figure 3H. Inset displays the ordered projection of Drosophila wing hair 
in proximal-distal axis. Arrow and hammer indicate the activating and repressing roles of kinase Fj against Ft and Ds, respectively. (B) Asymmetry of 
Ft-Ds heterodimer is disrupted when phosphorylation of Ft by Fj kinase is lost in Fj (Ft) mutant. This would result in swirling of wing hairs (inset). (C) 
Loss of kinase activity of Fj against Ds in Fj (Ds) mutant, however, does not influence asymmetric localization of Ft-Ds heterodimer. Orientations of wing 
hairs now remain intact (inset) like that of wild type wing (A). 
doi:1 0.1 371 /journal.pone.0097641 .g005 



-DsU 



(8) 



mechanisms - sub-cellular polarity of Ds" and that of Ffl- to be 
implicitly equal. However, in order to investigate whether they 
contribute equally in generating the sub-cellular asymmetry of Ft- 
Ds heterodimer, we introduced a weight factor /.lin the equation as 



The term F^(Z>^ +1 - 



represents the difference be- 
tween the distal and proximal sub-cellular concentrations of the 
heterodimer formed by Ffl of the l cell (with Ds" of(i+l)"' cell and 
with Ds" of cell, respectively). Therefore, the term 

Ft? (Ds" + , — Ds" j ) denotes sub-cellular polarity of Ff t or 
phosphorylated Ft in the i cell. By similar reasoning, the other 
term Ds"(Ft" i+l —F^_^, denotes sub-cellular polarity of Ds"'m i' h 
cell. Hence, it is evident that overall sub-cellular asymmetry of 
heterodimer in i"' cell, ACf 9 , are contributed by sub-cellular 
polarity of phosphorylated Ft, Ff t (Ds" + , — Dsf_ i ) , and that of 
unphosphorylated Ds, Dsf(Ff i+l -Ff^). 

While sub-cellular asymmetry of the Ft-Ds heterodimer has 
been shown experimentally [18,19], previous experiments did not 
resolve if sub-cellular polarity of Ds" and that of FtP contribute 
equally to sub-cellular asymmetry of the heterodimer. We thus 
used in silico mutations to investigate this crucial question. In all 
simulations till now, we have considered contributions of these two 



ACf = 



■^0+^(^+1-^-1)] ( 9 ) 



When contribution of sub-cellular polarity of Ds" in determin- 
ing overall sub-cellular polarity of the cell was taken to be much 
higher (ji> > 1 in the in silico mutations) than that of Ffl , we found 
sub-cellular asymmetry of heterodimers to be much more 
pronounced (Figure 6B) as compared to a condition when their 
contributions were considered to be equal - that is ft = 1 
(Figure 6A). This suggests that sub-cellular polarity of Ds" 
amplifies sub-cellular asymmetry of Ft-Ds heterodimer. On the 
contrary, when contribution of sub-cellular polarity of Ft was 
taken to be much higher, or that of sub-cellular polarity of Ds" was 
taken to be much smaller (jt< < 1), asymmetric localization of 
heterodimer was lost, or even reversed in some distal most cells 
(Figure 6C). This reinforces the interpretation that sub-cellular 
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polarity or asymmetric localization of Ds" largely contributes to 
sub-cellular distal enrichment of the Ft-Ds heterodimer. 

As explained above, Ft-Ds heterodimer at the distal edge of i h 
cell can be formed in two alternative orientations. However, it has 
not been resolved experimentally which of these two orientations 
(Ds"Ff i+x or FtfDs? +l ) is the preferred orientation of the Ft-Ds 
heterodimers at the all cell edges within a tissue. It is likely that 
these two alternative orientations of the heterodimer are linked to 
sub-cellular polarity of phosphorylated Ft and also that of 
unphosphorylated Ds. As mentioned above, sub-cellular polarity 
of unphosphorylated Ds, or Ds" in a cell is given by 
Ds"(Ff j+l —Ft p i X Y Similarly, sub-cellular polarity of FtP'm a cell 
is given by Ff t (Ds u i+ ; — Ds" , J . We found that overall sub-cellular 
asymmetry of the Ft-Ds heterodimer is amplified by sub-cellular 
polarity of unphosphorylated Ds (Figure 6B); therefore, the 
preferred orientation of the heterodimer at all cell edges in the 
tissue would be the one that corresponds to sub-cellular polarity of 
unphosphorylated Ds, i.e. Ds"Ff j+x and not to that of phosphor- 
ylated Ft, i.e. Ff i Ds" +l . Thus, our results indicate that heterodi- 
mer formed at the distal edge of a cell is between its unpho- 
sphorylated Ds and phosphorylated Ft of its distal neighbor [Ds" of 
z'Vll with Ff of (i+lf cell, or DsfFf i+l ). This presents an 
interesting possibility, namely, that within a cell, unphosphorylated 
Ds is asymmetrically localized towards its distal edge, while 



phosphorylated Ft is localized towards its proximal edge. 
Immunochemical detection using an antibody that discriminates 
between phosphorylated and unphosphorylated forms of Ft and 
Ds can help validate this prediction. 

In a recent phenomenological model of PCP [31], subcellular 
asymmetry of Ft-Ds heterodimer and a consequent tissue-scale 
collective cell polarization were achieved without considering the 
Fj gradient. Rather, Ft-Ds heterodimer asymmetry at the cell 
edges was obtained by assuming mutual inhibition between the 
two alternative orientations of the heterodimers while promoting 
their respective orientations. In other words, Ft p i Ds" +l promotes 
formation of more of its kind while inhibiting formation of 
DtffFf i+l ; likewise, Ds"Ft^ +i does the converse [31]. While this 
remains an attractive possibility, available experimental data do 
not offer support for such mutual inhibitory and self-promoting 
roles of Ft-Ds heterodimers. In contrast, our model does not 
assume such a mutually inhibitory and self-promoting role of Ft- 
Ds heterodimers and yet achieves their sub-cellular asymmetry at 
the cell edges by factoring Fj gradient (Figure 4A, C). A resolution 
of the question as to which one of these two models approximates 
an actual cellular event would demand experimental investigation. 
It may be further noted that this model proposed that all cell edges 
in a tissue would prefer any one out of two orientations of the 
heterodimers (FtfDs" +i and Ds"Ff i+l ) without a specific prefer- 
ence of one over the other [31]. Given that Fj kinase modulates the 
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Figure 6. Sub-cellular asymmetry of £>.s"enriches asymmetry of Ft-Ds heterodimer. (A) Equal contributions of sub-cellular asymmetries of 
Ft and Ds" result in asymmetric localization of Ft-Ds heterodimer (u = 1). Inset shows the enlarged view with inclined tops of the trapezoids. Data 
shown here is same as that shown in Figure 3H. (B) Higher contribution of sub-cellular asymmetry of Ds" augments the asymmetric localization of Ft- 
Ds heterodimer (steeper top of the trapezoids) in cells (u= 10). (C) Higher contribution of Ff sub-cellular asymmetry, however, results in diminished 
or even reversed asymmetry of Ft-Ds heterodimer in cells (|i = 0.1). Inset shows the enlarged view with flatter tops of the trapezoids indicating weak 
asymmetry of Ft-Ds heterodimer. 
doi:1 0.1 371 /journal.pone.0097641 .g006 



PLOS ONE | www.plosone.org 



8 



May 2014 | Volume 9 | Issue 5 | e97641 



Sub-Cellular Asymmetry of Ft-Ds Heterodimer for Planar Cell Polarity 



affinities of Ds and Ft for each other, a consideration of Fj gradient 
made in our model shows that Ds"Ft p +] , and not Ft?Ds" +l , is the 
preferred orientation of Ft-Ds heterodimer at the cell edges 
(Figure 6B, C). 

Caveats and Concluding Remarks 

A caveat of our model is its unidimensionality: interactions 
within the Ft-Ds-Fj system have been captured in only one 
dimension within the plane of tissue. Adopting this model in two 
dimensions presents challenges in terms of defining the direction of 
polarity. Further, our analysis explains sub-cellular asymmetry of 
Ft-Ds heterodimer; deciphering interactions of the global regula- 
tors (Ft-Ds-Fj system) with the core proteins (Fz, Vang, Pk, Dsh) 
remains outside the scope of our model. These caveats notwith- 
standing, our analysis shows that a mathematical model of PCP at 
an appropriate level of abstraction can provide novel and non- 
intuitive insights about the dynamics of the system [35] and 
testable predictions. For instance, our simulations predict a tissue- 
level spatial gradient of the Ft-Ds heterodimer concentration 
(Figure 3H). We also predict that sub-cellular polarity of Ds 
amplifies sub-cellular asymmetry of Ft-Ds heterodimer (Figure 6B, 
6C), thereby revealing a fundamental feature of a planar polarized 
cell. Similarly, we predict that phosphorylation of Ft by Fj is 
essential for establishing PCP, and suggest specific mutational 
studies to validate the same (Figure 5B, 5C). 

In conclusion, our mathematical model reveals that a set of 
known interactions within the Ft-Ds-Fj system, namely, gradient of 
Fj expression, phosphorylation of Ft by Fj kinase, and sub-cellular 
Ds localization generate an asymmetric enrichment of the Ft-Ds 
heterodimer in distal cell membrane. This could in turn provide 
PCP cue across the tissue. Our results are in agreement with 
previous experimental findings which reveal that Fj gradient, but 
not Ds gradient, is required for establishing PCP [16,33], and that 
PCP is more resistant to overexpression of Ds or Ft in the distal 
region of the wing as compared to the proximal [16,32]. It may be 
further noted that global regulators and core proteins display 
tissue-specific interactions. In Drosophila wing and eye, for instance, 
global regulators acts upstream of the core module [6] - the latter 
in turn act on tissue-specific effectors for PCP. In the abdomen, on 
the other hand, Ft-Ds-Fj system is proposed to signal directly to 
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